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ABSTRACT 

Gravitational microlensing surveys target very dense stellar fields in the local group. As a 
consequence the microlensed source stars are often blended with nearby unresolved stars. The 
presence of 'blending' is a cause of major uncertainty when determining the lensing proper- 
ties of events towards the Galactic centre. After demonstrating empirical cases of blending 
we utilize Monte Carlo simulations to probe the effects of blending. We generate artificial 
microlensing events using an HS T luminosity function convolved to typical ground-based 
seeing, adopting a range of values for the stellar density and seeing. Microlensing light curves 
are generated using typical sampling and errors from the second phase of the Optical Gravi- 
tational Lensing Experiment. We find that a significant fraction of bright events are blended, 
contrary to the oft-quoted assumption that bright events should be free from blending. We 
probe the effect that this erroneous assumption has on both the observed event timescale 
distribution and the optical depth, using realistic detection criteria relevent to the different 
surveys. Importantly, under this assumption the latter quantity appears to be reasonably unaf- 
fected across our adopted values for seeing and density. The timescale distribution is however 
biased towards smaller values, even for the least dense fields. The dominant source of blend- 
ing is from lensing of faint source stars, rather than lensing of bright source stars blended 
with nearby fainter stars. We also explore other issues, such as the centroid motion of blended 
events and the phenomena of 'negative' blending. Furthermore, we breifly note that blending 
can affect the determination of the centre of the red clump giant region from an observed 
luminosity function. This has implications for a variety of studies, for example mapping ex- 
tinction towards the bulge and attempts to constrain the parameters of the Galactic bar through 
red clump giant number counts. We conclude that blending will be of crucial importance for 
future microlensing experiments if they wish to determine the optical depth to within 10 per 
cent or better 
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1 INTRODUCTION 

Gravitational microlensing is maturing into an important astrophys- 
ical technique with diverse applications to Galactic astronomy, such 
as probing the dark matter content of the inner Galaxy (see, e.g., the 
following review articles: Paczynski 1996; Mao 1999; Evans 2003). 
Thousands of microlensing events have been discovered. The vast 
majority of these are towards the Galactic centre and many were 



identified in real-time, for example by the OGLlfl (Udalski 2004) 
or (Bond et al. 2001) alert systems. The microlensing prob- 

ability (known as the optical depth, r) towards the Galactic centre 
probes the mass distribution along the line of sight. The earliest de- 
terminations yield optical depths (Udalski et al. 1994a; Alcock et 
al. 1997a; Alcock et al. 2000) that are significantly higher than the 
theoretical predictions (e.g., Zhao & Mao 1996; Binney, Bissantz 
& Gerhard 2000; Evans & Belokurov 2002; Han & Gould 2003). 
More recent determinations yield both lower values (Popowski et 
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al. 2005; Hamadache et al. 2006; Sumi et al. 2006) and higher val- 
ues (Sumi et al. 2003), although these determinations are based on 
relatively small samples of microlensing events. It is important to 
note that there appears to be a clear distinction between the mea- 
sured values of t for the two commonly-used techniques: higher 
values of t are found for determinations carried out using all stars 
(e.g. Alcock et al. 2000; Sumi et al. 2003), whereas lower values 
are found when (as advocated by Gould [1995]) only bright stars 
are used in the analysis (e.g. Popowski et al. 2005; Hamadache et 
al. 2006; Sumi et al. 2006). 

It was realised quite early on (e.g. Udalski et al. 1994a; Al- 
cock et al. 1997b) that blending is a major uncertainty in the de- 
termination of T. Blending occurs naturally because microlensing 
surveys are conducted in crowded stellar fields and, with typical 
ground-based seeing, other stars can blend into the seeing disk of 
the lensed star (see, for example Han 1999 and references therein). 
This affects the number of potential lensed sources and also intro- 
duces uncertainties into the fitted event parameters (e.g. Wozniak 
& Paczyhski 1997; Han 1999). Importantly, it was proposed that 
the aforementioned discrepancy between the t measurements from 
all stars compared to bright stars could be explained by blending. 

Clearly the ideal way to understand this blending issue is 
with high resolution HubbleS paceTelescope(HST) images (Han 
1997), which can be used to resolve any nearby blends that may 
be present. This technique was adopted by Alcock et al. (2001a) 
in their analysis of a set of microlensing events towards the Large 
Magellanic Cloud. However, in general this method is limited due 
to the restrictions on the availability of HS T observing time. There- 
fore, in the absence of high resolution images for each event, the 
next best approach is to undertake Monte Carlo simulations at the 
pixel level. Mock images can be generated based on deep HS T lu- 
minosity functions of the Galactic bulge. Unfortunately, since such 
luminosity functions are currently only available for a very limited 
number of lines of sight (e.g. Holtzman et al. 1998), one must ex- 
trapolate their behaviour for the various bulge fields observed by 
microlensing collaborations. Artificial microlensing stars can be 
injected into mock images and then convolved into ground-based 
seeing. It is then possible to investigate the efficiency of recover- 
ing microlensing events, i.e. the detection efficiency (see, for ex- 
ample, Alcock et al. 2000). In order to simplify the analysis many 
microlening studies concentrate on bright stars, working under the 
assumption that bright stars suffer negligible blending. However, 
the reliability of this assumption has recently been called into ques- 
tion (see Section|2ll. As a result, proper consideration must be given 
to blending, even when one considers microlensing of bright stars. 

Monte Carlo simulations of blending have already been car- 
ried out by various groups, mostly concentrating on the effect to 
the recovered microlensing optical depth. Recent studies include 
Popowski et al. (2005) and Hamadache et al. (2006), both of which 
argue that the recovery of the optical depth should not be signifi- 
cantly biased by the presence of blending in bright events. The anal- 
ysis presented in this paper builds upon another such work (Sumi 
et al. 2005), which showed that a simulated sample of bright mi- 
crolensing events still contains many heavily blended events. We 
extend the work of Sumi et al. (2005) by generalising the analysis 
to fields with varying stellar density under diff^erent seeing condi- 
tions. By doing this we aim to make broader conclusions that go 
beyond any experiment-specific analysis. 

The outline of the paper is as follows, in Section|2]we briefly 
discuss some of the observational evidence that exists to suggest 
that bright microlensing might not avoid the problem of blending. 
The remaining sections deal with Monte Carlo simulations that we 



have undertaken in order to investigate this phenomenon: Section|3] 
describes the method; Sections|4]&[5]present the results of our sim- 
ulations including the resulting distributions of event parameters; 
Section|6]investigates whether the assumption that bright events are 
unblended can bias the measured value for the optical depth; and 
Section|7]discusses the expected distributions of centroid motions. 
We conclude with Section[8l where we discuss the implications of 
our findings. 



2 OBSERVATIONAL EVIDENCE OF BLENDING IN 
BRIGHT EVENTS 

Although it is often assumed that in general bright events are not 
affected by blending, there is observational evidence to show that 
this is not always a safe assumption. In the following section we 
briefly discuss various different approaches that can be used to test 
this hypothesis for observed events using ground-based data. We 
provide a number of examples of bright events that exhibit signifi- 
cant blending. 

Throughout this section and the rest of the paper we charac- 
terise the blending using the following parameter, /s, which denotes 
the ratio of the source flux to the total baseline flux, i.e. 

p ^ source ^ 

Is = . (1) 

-•^baseline 

Therefore, /s = 1 for the case of no blending, while /s — > for 
heavily blended events. 

2.1 General model fitting 

Numerous catalogues of gravitational microlensing events have 
been published towards the Galactic bulge. Many of these con- 
tain bright events that model fitting has suggested are blended, such 
as the MACHO catalogue of Alcock et al. (2000). This catalogue 
contains a subset of 37 events that are described as 'classical lens- 
ing' from which one can determine the fitted blending parameter. 
Four of these events are both bright (V < 18) and heavily blended 
(/s < 0.2): namely 95-BLG-dl9, 97-BLG-dl3, 97-BLG-24, 97- 
BLG-37. All of these four events have reasonably well constrained 
values of the blending parameter. Popowski et al. (2005) also noted 
that based on light curve fitting they cannot exclude the possibility 
of significant blending for some events in their clump-giant sample. 

Further examples can be found in the catalogue of clump-giant 
EROS events published by Afonso et al. (2003), in which two of 
16 bright events were found to display clear blending signatures 
(EROS-BLG-31 and EROS-BLG-12). In contrast, a more recent 
analysis of the EROS data (Hamadache et al. 2006) finds only five 
of 120 clump giant events appear to exhibit strong blending, al- 
though they acknowledge that their paucity of blended events could 
be due to the limited photometric precision of the EROS experi- 
ment. 

In their analysis of the OGLE-II catalogue of bright events, 
Sumi et al. (2006) found a that blending was significant for a num- 
ber of their events. According to their best-fitting models, 38 per 
cent of these bright events were actually due to lensing of much 
fainter sources. 

2.2 Centroid motion 

Another way to assess the influence of blending is to investigate 
the motion of the light centroid during the event. The light cen- 
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troid is determined by the liglit from the lens, the lensed source, 
and/or blended stars along the line of sight. Therefore, during mi- 
crolensing, the light centroid must shift towards the lensed star as 
it brightens (e.g., Alard, Mao & Guibert 1995; Goldberg 1998). If 
the lens dominates the blend, the centroid shift is difficult to detect 
because the lens and the lensed source are aligned to within mil- 
liarcseconds (~ angular Einstein radius). However, if the blending 
is mostly due to other blended stars, then the centroid shift may be 
detectable, even in ground-based observations. 

Another similar approach is to measure the offset between the 
lensed source and the centroid of the light at baseline (Han 2000); 
this approach can be used when dealing with Difference Image 
Analysis (DIA), since in this case the location of the lensed source 
can be measured to high accuracy. It is also possible to try and 
locate the blend by removing the light of the lensed source using 
image subtraction (Gould & An 2002; Smith et al. 2002). 

These techniques have been applied in various works. The off- 
set between the lensed source and the baseline centroid has been 
routinely measured for various microlensing catalogues (e.g. Al- 
cock et al. 1999; Alcock et al. 2000; Wozniak et al. 2001). In many 
cases there are a significant fraction of events with offsets of 1 arc- 
sec or greater. 

The first such detections of centroid motions have been pre- 
sented in Alard, Mao & Guibert (1995) and Goldberg & Wozniak 
(1998). To further test this effect we have examined the centroid 
motion for a sample of red clump giant microlensing events from 
OGLE-II (Sumi et al. 2005). In order to constrain the offset between 
the lensed star and the blend, we fit the light curve and the centroid 
positions simultaneously. Fig. [T] shows the event with the most sig- 
nificant centroid shift (sc37-556534), which has a very bright base- 
line magnitude (/ = 15.9 mag). As expected, the centroid location 
is a strong function of the magnitude (magnification). The offset 
(in pixels, where 1 pixel corresponds to 0.417 arcsec) between the 
lensed star and the blend is hx a 0.61, Ay a -1.04. This centroid 
shift clearly demonstrates that blending can be important for bright 
microlensing events. 



2.3 Binary lens blending 

Blending can be inferred not only for single microlensing events, 
but also for binary microlensing events. In fact, the first (Udal- 
ski et al. 1994b) and second (Alard, Mao & Guibert 1995) bi- 
nary lens events both show significant blending. Detailed studies of 
binary lenses in OGLE-II (Jaroszynski 2002) and OGLE-III data 
(Jaroszynski et al. 2004) show convincingly that blending in binary 
lenses is widespread; they find that the fractions of bright events 
(observed / < 18) with fitted parameter /s < 0.5 are 3/5 and 4/7 
for OGLE-II and OGLE-III, respectively. The most spectacular ex- 
ample is sc5_6650 from OGLE-II, for which the lensed source is 
inferred to contribute only I per cent of the total light even though 
the composite is very bright with a baseline magnitude / = 16.18. 

Binary light curves are very diverse and, as a result, poorly 
sampled ones can often be fitted with multiple models (e.g., Mao 
& Di Stefano 1995; Albrow et al. 1999; Dominik 1999; Gaudi & 
Han 2004). This problem often renders the /s determination some- 
what uncertain. However, for binary events that undergo a caustic 
crossing, a limit on the blending can often be inferred without any 
detailed modelling. These binary events exhibit 'U'-shaped light 
curves as the lensed star enters and then exits from the caustics. 
The minimum magnification in the plateau must be equal to or ex- 
ceed 3 (Witt & Mao 1995). If the observed minimum magnification 



within this plateau, A^in, is below 3, then the fraction of light con- 
tributed by the source must satisfy the inequality, 

h < A™„ ^ 3. (2) 

For example, the OGLE-II event sc5_6650 has an observed A,ni„ « 
1.04 and hence we can infer that /s < 0.02, which is fully consistent 
with the blending parameter (/s = 0.01) obtained by Jarosznski 
(2002) from detailed light curve fitting. Kim (2004) presents a more 
comprehensive analysis of the limits that can be derived for binary 
events in the OGLE-II and OGLE-III databases; this work shows 
that 5/7 bright (observed / < 18) caustic crossing events have /s < 
0.5. 

It should be noted that there may be a slight bias in the blend- 
ing distribution derived from caustic crossing binary events since 
such events often undergo large amplifications, which increases 
the probability of observing lensing of faint source stars. However, 
even given this caveat it is clear that the analysis of binary events 
can provide a robust model-independent method for investigating 
blending. As has been shown above, it is evident that for binary 
lenses (as with single lenses) bright events can be heavily blended. 



3 CONSTRUCTION OF A MOCK CATALOGUE OF 
BLENDED EVENTS 

In the remainder of this paper we undertake Monte Carlo simu- 
lations in order to investigate the effect of blending for simulated 
bulge fields with varying seeing and densities. In this section we 
discuss the construction of our catalogue of mock microlensing 
events. 

We first construct artificial images for 9 simulated fields, 
adopting three different values for the density of stars and three 
different values for the level of seeing. The densities of our fields 
are chosen relative to the OGLE-II field sc3, which is centered on 
/ = 0.11°, b = -1.93° and has an observed surface density of 151 
stars per arcmin^ down to a magnitude of / = 17. We choose densi- 
ties of 0.5, 1 and 1.5 times the density of field sc3. The three values 
of seeing are 2.1, 1.05, 0.7 arcsec. Throughout this paper we desig- 
nate the field with median seeing and density (i.e. density of field 
sc3 and seeing of 1.05 arcsec) as our reference field. The details for 
each of our fields are summarised in Table[T] Similar to the OGLE- 
II experiment, each field has dimensions 14x57 arcmin- and has 
pixel size 0.417 arcsec pixeL'. Field sc3 was chosen since it is a 
very dense stellar field close to the Galactic centre where blend- 
ing should be most significant. Further details about the OGLE-II 
experiment can be found in Udalski et al. (2000). 

It is helpful to see how the adopted characteristics for our sim- 
ulated fields compare to the important microlensing bulge surveys. 
The pixel size and typical seeing values for five of the major mi- 
crolensing experiments are: 

• OGLE-II: 0.42 arcsec pixel size and median seeing ~ 1.3 arc- 
sec (Sumi et al. 2006) 

• OGLE-III: 0.26 arcsec pixel size and median seeing ~ 1.3 arc- 
sec 

• MACHO: 0.63 arcsec pixel size and median seeing ~ 2. 1 arc- 
sec (Popowski et al. 2005) 

• EROS-II: 0.6 arcsec pixel size and median seeing ~ 2 arcsec 
(Palanque-Delabrouille et al. 1998); 

• MOA: 0.81 arcsec pixel size and median seeing ~ 2.5 arcsec 
(Bond et al. 2001). 
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Figure 1. An example of an observed case of centroid motion. The top panel shows the light curve for this event (sc37-556534). The bottom panel shows the 
corresponding centroid shift for ,y and v in pixels (1 pixel corresponds to 0.417 arcsec) as a function of /-band magnitude. The solid line indicates the best-fit 
model. For clarity, the y centroid position has been shifted downward by 1 pixel. 



It can be seen that the range of seeing and density for our simulated 
fields covers practically all bulge surveys, both current and pre- 
vious. Although our adopted pixel size (corresponding to OGLE- 
II) is smaller than the non-OGLE experiments listed here, this 
should not be a dominant effect. It is also worth noting that our 
choice of medium density, which corresponds to one of the densest 
OGLE-II fields, is comparable to the densest fields in both MA- 
CHO (Popowski et al. 2005) and EROS (Hamadache et al. 2006) 
optical depth studies. 

Our artificial fields are created following the prescription of 
Sumi et al. (2006); full details of the procedure can be found in sec- 
tion 4. 1 of their paper. We populate our field by selecting stars from 
the HST /-band luminosity function of Holtzman et al. (1998). 
Since the OGLE field sc3 is not coincident with the HST field 
from Holtzman et al. (1998), the HST luminosity function must 
be shifted so as to matches the observed number density of bright 
stars from the OGLE sc3 field. By combining the OGLE and HS T 
data we are also able to constrain the bright end of the luminos- 
ity function; this region cannot be well constrained using the HS T 
data alone as there are very few bright stars due to saturation and 
the small field of view. We account for the differential extinction 
across the field using the extinction maps of Sumi (2004). Note 
that our simulations are carried out solely with /-band data and do 
not incorporate any colour information. 

Using this luminosity function we populated our field with ap- 
proximately 10^ artificial stars down to a magnitude of / « 22 and 
then convolved the image to our required level of seeing. Our re- 



sulting mock fields are therefore fully synthetic (i.e. artificial stars 
are not injected into observed fields) and have realistic noise prop- 
erties, making them effectively indistinguishable from real images. 
Note that stars are placed randomly in the field and are not placed 
on a regularly spaced grid. 

Given this mock field we then applied the standard OGLE 
photometry pipeline that is based on DoPHOT (Schechter, Mateo, 
& Saha 1993) to obtain the 'observed' magnitude of the lO' arti- 
ficial stars. From this we can record the 'input' magnitude of the 
artificial star (/jn) and the resulting 'output' magnitude as measured 
by DoPHOT (/out). This 'output' magnitude corresponds to the ob- 
served baseline for a star. Ideally, we should find that Ia„i < /in for 
all stars, since the 'output' magnitude includes many blended stars. 
However, as will be shown below, this is not true in all cases (see 
Sectionl42t. 

Given these 'input' and 'output' magnitudes we can calculate 
the the blending fraction, /s, i.e. the fraction of the baseline flux 
contributed by the source, 

j.^ _ -^soui-ce _ -^source _ j^Q(;„u,-/i„)/2.5 p-j 

^baseline ^source ^blend 

However, the resulting blending distribution needs to be corrected 
for the detection efficiency, since heavily blended events (i.e. fs ^ 
1) are less likely to be detected because they require much greater 
intrinsic magnifications to produce an observed increase in magni- 
tude. 

To simulate this effect, we generated a mock catalogue of stan- 
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dard microlensing light curves using the sampling and photomet- 
ric properties of the OGLE-II experiment (Udalski et al. 2000): 
namely, /-band observations taken approximately once every few 
nights; bulge season typically lasting from mid-February until the 
end of October; limiting magnitude and saturation are approxi- 
mately / =s 20 and / ~ 11.5, respectively. We choose to generate 
light curves over a total baseline of three years and assume that each 
star has an equal probability of being lensed (i.e. we assume all stars 
belong to the bulge). The generation of these mock catalogues fol- 
lows the prescription given in Smith et al. (2005) although for the 
purposes of this paper parallax signatures (Gould 1992) have been 
neglected. Although it is well known that symmetric parallax sig- 
natures can be confused with blending (Smith, Mao & Paczyhski 
2003), it would be computationally too demanding to generate and 
fit parallax signatures. In any case, such symmetric parallax events 
should not be common. 

The event timescale is drawn from the model described in 
Smith et al. (2005) and the impact parameter (in units of the Ein- 
stein radius) is chosen uniformly between and 1.5. Note that un- 
less otherwise mentioned, our analysis is based on events with fit- 
ted impact parameter less than 1; we generate events with larger 
impact parameters because when events are fitted with a microlens- 
ing model (see Section [STTTT) the recovered impact parameter can 
sometimes be underestimated. Photometric errors are assumed to 
be Gaussian (e.g. Wozniak 2000). 

Once this mock catalogue has been produced we apply the se- 
lection criteria of Sumi et al. (2005) to simulate the detection effi- 
ciency. Essentially, these criteria test for the presence of a constant 
baseline with one distinct brightening episode. 

We divide our mock catalogue into two samples: a full sample 
of all events with baseline magnitude /om < 19 and a subsample of 
bright events with < 17. This latter subsample of events corre- 
sponds to what one would expect to see for a sample of red clump 
giant microlensing events, for example. We only consider events 
with /out < 19 since fainter objects are too close to the limiting 
magnitude of the OGLE-II experiment. These cuts are illustrated 
in Fig. [2] which shows the observed colour-magnitude diagram for 
the OGLE-II sc3 field. 



4 GENERAL SHAPE OF BLENDING DISTRIBUTION 

4.1 General description of blending distribution 

The resulting blending distributions for 5 of our simulated fields 
are shown in Fig. |3] In this figure we show both the full samples 
of all events (/„„, < 19) and the bright subsamples of events with 
baseline magnitude brighter than 1„m < 17 (see Section[3j. 

For the full sample of all events, the distributions appear rel- 
atively flat between /s = and /s = 1. The most striking feature 
is the long tail of events with /s > 1, particularly for fields with 
bad seeing. This feature corresponds to sources that are brighter 
than the measured baseline (i.e. /om). These so-called 'negatively 
blended' events are discussed in detail in Section |4!2l As expected, 
there is a significant proportion of heavily blended events; for ex- 
ample, for our reference field (medium density and seeing), 57 per 
cent have /s < 0.5 and 16 per cent have /s < 0.1. For our dif- 
ferent fields the blending distribution follows the anticipated trend, 
with the better seeing and/or less dense stellar fields exhibiting less 
blending (see Table[T]l. 

The shape of the blending distribution is quite different for the 
bright subsample of events. There is a clear 'u' -shaped distribution 



with few events in the region 0.2 < /s < 0.8. It has been assumed 
that bright events may be free from the problem of blending, since 
any blended star would be too faint to contribute substantially to 
the total flux. However, it is clear from this figure that a significant 
proportion of events have /s < 0.1, namely 27 per cent for our 
reference field. The behaviour for the other fields is very similar, 
but with fewer heavily blended events for fields with better seeing 
and/or lower stellar density (see Table[Tll. 

We find that in almost all cases heavily blended events occur 
in our bright sample when a faint source close to a bright star is 
lensed. For example, in our reference field we find that for a sam- 
ple of bright (total baseline magnitude /„„, < 17), blended events 
(/s < 0.2), over 99 per cent have source magnitudes fainter than 
17th magnitude. There are two competing effects that will deter- 
mine how frequently such events occur. First, it depends on the 
luminosity function, i.e. how the number counts rise as the magni- 
tude becomes fainter. If the luminosity function is steep (as is the 
case for the bright end of the luminosity function), then there are 
numerous faint stars, which enhances the number of lensed cases 
with strong blending. Second, in order for strongly blended events 
to be observable, the intrinsic magnification must be much higher 
than the nominal threshold (A > 3/ VS). Given the shape of this /s 
distribution, it appears that the assumption that bright events are un- 
blended is not a sound one. We deal with the effects of this assump- 
tion in later sections (see Section[5]and[6ll. As with the full sample 
of all events there is a fraction of events with negative blending, al- 
though for this magnitude range this fraction is smaller (see Section 

In Fig. |4] and Table |2] we show how the /s distribution varies 
as a function of the baseline magnitude (i.e. /o^) for our reference 
field. Although the overall shape of the distribution shows a signifi- 
cant trend across the range of magnitudes, there is no clear trend for 
the percentage of heavily blended events with /s < 0.2. This fact 
is interesting in that it is clearly in contradiction with the argument 
that bright events are free from significant blending. 

4.2 Negative blending 

Positive blending is obviously induced when the lensed source 
is blended with other sources that are too close to resolve under 
ground-based seeing conditions. The case for negative blending is 
less straightforward. One of the first papers to discuss the issue of 
negative blended fluxes was Park et al. (2004), which concerned the 
microlensing event MOA-2003-BLG-37. For this event they obtain 
a best-fit value of /s ~ 1.06. Further events have been identified by 
various authors (e.g. Jiang et al. 2004; Poindexter et al. 2005; Sumi 
et al. 2006). 

Fig. [3] shows the extent of the issue for our simulated fields 
(see also Table [T). Curiously, there are no obvious trends evident 
in either seeing or density. The only clear effect is that for lowest 
quality seeing the problem becomes significantly worse, although 
this could be due to increased photometric noise which is present 
in fields with such bad seeing. There is also no clear trend between 
the samples of all events and bright events. 

We examined a number of images from our simulations that 
show negative blending. In many cases we find that it arises when 
we have another star close to the lensed star, and the DoPHOT pho- 
tometry program incorrectly partitioned part of the lensed star flux 
to the nearby blend, yielding a negative blending for the lensed 
source. Another possibility is a 'hole' in the mottled distribution 
of faint stars that constitutes the background flux. This arises be- 
cause the DoPHOT program assumes a constant background, while 
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Figure 2. The observed OGLE-II colour-magnitude diagram for field sc3. Note that this is the observed field and not a simulated field, since we do not include 
V-band data in our simulations. The horizontal lines show the cuts for all events {laM < 19) and bright events (/om < 17)- The latter sample coiTesponds to a 
brightness cut similar to that employed for red clump giant microlensing events; for this figure the red clump region is centred around V - I ^ 2.4, / ^ 16. 
In reality, a red clump sample would be selected using a cut in colour to remove foreground blue disk stars, but this is not necessary in our case since our 
simulation contains no colour information and is set up so that all stars belong to the bulge. 



in reality the background is contributed by unresolved faint stars 
which have fluctuations. If the background close to the star is lower 
than the average background (a 'hole'), then the DoPHOT program 
will over-subtract the background, yielding a negative blending for 
the lensed source. For MOA-2003-BLG-37, the best-fit deblended 
intrinsic source magnitude Ig a 14.4 mag, in order to produce 
/s ~ 1.06, Park et al. (2004) required a 'hole' equivalent to the 
omission of a /q a 17 mag turn-off star. 

To conclude, we are unable to make any definitive statements 
about the nature of negative blending. We have attempted to show 
the empirical effects by carrying out analyses that are as close to 
real surveys as possible. It seems that the causes for negative blend- 
ing are a combination of those discussed above, namely software 
issues from deblending and holes in the background, with the addi- 
tional complication of simple statistical noise. 



5 INVESTIGATION OF BLENDING-RELATED BIASES 
IN THE RECOVERY OF EVENT PARAMETERS 

5.1 Fitting of the mock light curves and detection efficiency 

5.1.1 Light cun>e fitting 

To investigate the best-fit event parameters we fit each mock light 
curve with the standard 5-parameter blended microlensing model: 



/(f) = /base -2.5 log [/s(A(f)-l)+l], 

where the magnification, A(t), is given by. 



A(t) = 



u(tY + 2 
u(t) ^ju{ty + 4 



[((f) = 



(4) 



(5) 



Here uq is the impact parameter in units of the Einstein radius, fo 
is the time of the closest approach (i.e. maximum magnification), 
and % the event timescale. Our parameter % corresponds to the 
Einstein radius crossing time; it should not be confused with the 
Einstein diameter crossing t = It^, which is sometimes used by the 
MACHO collaboration. 

The best fit was found by minimizing the using the Minuit 
package from the CERN Program Library|f| A full description of 
the Minuit package can be found in the 'Minuit Reference Manual' 
(James 1994), which is also available online. 

We reject all events for which the timescales are degenerate 
(i.e. those events for which errors cannot be computed by Minuit or 
events with errors on t^ greater than 30 per cent), since no meaning- 
ful parameter values can be extracted from such light curves. This 
mostly removes low signal-to-noise ratio (S/N) events, although a 



http://wwwasdoc.web.cem.ch/wwwasdoc/minuit/ 
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Figure 3. The input (i.e. not fitted) blending distributions for recovered microlensing events in our simulated fields. The upper panel shows the distribution 
for a bright subset of events (/oui < IV), such as red clump giants, and the lower panel shows the distribution for all events (i.e. /om < 19). The left panel shows 
fields with differing stellar density (dotted/solid/dashed correspond to high/medium/low density) while the right panel shows fields with differing values for 
the seeing (dotted/solid/dashed corresponds to low/medium/high quality seeing). The vertical line at /s = 1 denotes zero blending, which means that all events 
to the right of this line exhibit negative blending (see Section|4j2). The small-scale fluctuations in these distributions are due to small number statistics and are 
not physical. 



fraction of these are heavily blended events suffering from the well- 
known degeneracy between /s, uo and (Wozniak & Paczyiiski 
1997). 

Note that it is important that we only deal with events with 
well-constrained parameters since when we calculate the optical 
depth for the blended fits (Section[6j we select our bright samples 
using the fitted source magnitude, which obviously depends on /g. 

5.1.2 Notation convention 

In the remainder of this work we adopt the following notation con- 
ventions. We use the subscript 'in' to denote an 'input' property, i.e. 
the true value of a parameter. For events that have been fitted with 
a model incorporating blending (i.e. a 5-parameter fit) we use the 
subscript '5p'. Equivalently, for a fit with a model that incorporates 
no blending (i.e. a 4-parameter fit with /s = 1) we use the subscript 
'4p'. The subscript 'out' is sometimes used without an accompany- 
ing '4p' or '5p' to denote the fitted parameter in general, i.e. when 
referring to both 4- and 5-parameter fits. 

For the optical depth Tin corresponds to the input r, as given 
below by equation l|6j, and Tout corresponds to the observed value 
obtained using the best-fitting t^ parameters for the sample of mi- 
crolensing light curves (equation [7] below). 

We also retain the convention for magnitudes described above. 



namely that /;„ refers to the magnitude of a star on the 'input' im- 
age (i.e. prior to convolving to ground-based seeing) and /o„t refers 
to the magnitude of a star on the 'output' image (i.e. the result- 
ing magnitude of a star from DoPHOT after the image has been 
convolved to ground-based seeing). The magnitude /s denotes the 
magnitude of the source star, i.e. the star that has been microlensed. 

5.1.3 Detection efficiency 

The detection efficiency is required when one wishes to compare 
timescale distributions or calculate optical depths (see Section|6]l. 
It is is simply the proportion of microlensing events that are re- 
covered as a function of the event timescale, although it can be 
extended so that it becomes a function of other parameters such as 
the peak magnification (e.g. Popowski et al. 2005). We only require 
the efficiency for bright events with I^^, < 17 mag, since we will 
not deal with the fainter events in such a rigorous way. 

We calculate the detection efliiciency separately for our 
blended and unblended fits. For our unblended fits we are working 
under the (not necessarily correct) assumption that events brighter 
than /out = 17 mag are unblended. Therefore for each of our sim- 
ulated fields we generate a series of mock unblended events using 
the same method as described in Section[3] where the baseline mag- 
nitude of the source is selected from our standard luminosity func- 
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Figure 4. The input (i.e. not fitted) blending distribution for recovered microlensing events for six different output magnitudes, /om- This figure shows the 
distribution for the simulated field with medium seeing and density. The small-scale fluctuations in these distributions are due to small number statistics and 
are not physical. 



tion. We then apply our 4-parameter fitting routine (Section [5.1.1t 
and calculate the fraction of events that pass our detection criteria 
as a function of event timescale. In this case the detection criteria 
are based on the 4-parameter fits, e.g. ;(o.4p < 1. 

The detection efficiency for blended events is very similar 
to the unblended case, except for these we also incorporate the 
blended component of the baseline flux into our event generation. 
This is done using the simulations described in Section[3] We then 
fit our events with the 5-parameter fit and calculate the proportion 
of events that pass our detection criteria. Note that the efficiency is 
much reduced owing to the fact that many of these blended events 
are degenerate and thus do not pass our selection criteria. It is also 
important to note that for the blended fits we classify bright events 
as those with source magnitude brighter than /s = 17 mag, which 
means that events with < 17 can be excluded if the are suffi- 
ciently blended. 



5.2 The effect of blending for all events 

The distributions of the fitted and input event parameters are shown 
in Fig. [5] for our reference field (medium stellar density and see- 
ing). The bottom panels show the distributions for all events, i.e. 
/out < 19. For the blending parameter /s, the fitted distribution 
matches the input distribution reasonably well except for a slight 
under-prediction of events with /s sp just less than one and a slight 
over-prediction for the number of negatively blended events with 



/s,5p > 1. Note that this /s.sp distribution differs from that shown 
in Fig. [3] since this does not include degenerate events from which 
one cannot gain a reliable estimate on the event parameters. This 
leads to a reduction in the number of events with small /s.sp ow- 
ing to the well known degeneracy between /s, uo and (Wozniak 
& Paczyfiski 1997). In addition, such heavily blended events often 
have lower signal-to-noise owing to the fact that these faint source 
stars require larger amplifications to rise above the detection thresh- 
old. 



The bottom-right panel shows the fitted timescale distribution 
for all events {la^t < 19). The input timescale distribution (%_]„) 
for the events that pass the microlensing selection criteria almost 
precisely matches the distribution of fitted timescales (?E,5p)- The 
difference in (%) between these two distributions is less than one 
per cent. Note that we do not compare the behaviour of the dis- 
tribution with an actual observed distribution, such as in Alcock et 
al. (2000), since the nature of the timescale distribution is subject 
to various uncertainties which we do not wish to consider here. The 
purpose of this work is solely to show the difference between the 
fitted and input distributions. 

Although in Fig.[5]we only show the distributions for our ref- 
erence field, the other fields from our simulations exhibit qualita- 
tively similar behaviour. 
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Figure 5. Input and fitted event parameter distributions for the blending fraction (left) the event timescale (right). This figure corresponds to the simulated field 
with medium stellai' density and seeing; the other fields show qualitatively similar behaviour The solid line represents the input event parameter distribution, 
the long dashed line represents the 5-parameter fit (i.e. unconstrained blending) and the dotted fine represents the 4-parameter fit (i.e. fs = 1). The top panels 
show the distributions for a subset of bright events (/om < 17), such as red clump giants, and the bottom panels show the distributions for all events (i.e. 
/out < 19). The fitted timescale distributions in the upper panel are corrected for the detection efficiency and the solid line denotes the underlying timescale 
distribution; the lower panel simply shows the distributions of input and fitted timescales for the detected events (note that since the two distributions in this 
panel are almost identical, the two lines are coincident). 



5.3 The effect of blending for bright events 

The distributions of fitted event parameters for the bright sample 
of events from our reference field are shown in the upper panels of 
Fig-H] As was seen in Fig.|4l the distribution of blending parameter 
/s,in shows a 'u' -shaped distribution for bright events. However, the 
fitted value of the /s sp parameter shows a markedly smoother dis- 
tribution, where an under-prediction in the number of events with 
/s,5p ~ 1 results in an over-prediction of events with /s_5p ~ 0.5 and 
/s,5p > 1- 

Despite this problem, the distribution of f^.sp matches the in- 
put distribution extremely well. In the upper-right panel we show 
the fE,5p distribution corrected for the microlensing detection effi- 
ciency (see Section |5.1.3l >. along with t^^m distribution that we used 
to generate our events. The two distributions are in almost perfect 
agreement, with the fitted (?e) being recovered to within one per 
cent. 

It is interesting to see what happens to the fitted distribution 
when we impose the assumption that bright events are unblended. 
To investigate this we fit our sample of blended mock light curves 
with a model that enforces no blending (i.e. /s = 1), again correct- 
ing for the detection efficiency (Section [5. 1.3| l. The upper-middle 
panel of Fig. [5] shows the effect of this assumption. There is clearly 



a significant tail of short duration events in the distribution of fE,4p 
that is not present in the input distribution, which is due to the 
presence of blended events; when blended events are fitted with 
an unblended model, the resulting timescales are underestimated. 
However, the fit that incorporates blending does not suffer from this 
problem. Because of this problem, the (rg ^p) is under-estimated by 
23 per cent for this field. In the next section we discuss the impli- 
cations of this assumption on the optical depth. 

The qualitative behaviour for the other simulated fields is sim- 
ilar to that shown for our reference field. However, the effect of the 
assumption that bright events are unblended varies between fields, 
due to the fact that the level of blending is dependent on the stel- 
lar density and seeing. In Figure |6] we show the distribution of 
the bright sample of events for a selection of fields from our sim- 
ulations. This figure illustrates that the fields with a greater level 
of blending (i.e. those with greater density or worse seeing) have 
more problems reproducing the 1e distribution under the assump- 
tions that bright events are unblended. The corresponding depen- 
dence of {fE_4p) with density and seeing are tabulated in Table [T] 
and plotted in Figure|7] Note that while the value of (fE.oui>/<'E,in) 
is not dependent on the seeing or density for the blended fits (for 
all fields /e is recovered to within 2 per cent), there is a significant 
dependence for the unblended fits. 
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Figure 6. Input and fitted event timescale distributions for various values of density (bottom) and seeing (top). Note that these distributions are corrected for 
detection efliciency. The thin curves denote the input timescale distribution and the thick curves denote the fitted timescale distributions for the 5-parameter 
(i.e. unconstrained blending; dark) and 4-parameter (i.e. /s = 1; light) fits. In the bottom panels the dotted/solid/dashed curves correspond to high/medium/low 
stellar density, while in the top panels the dotted/solid/dashed curves correspond to low/medium/high quality seeing. All curves are only for bright events 
(/out < 17). 



One issue that concerned Hamadache et al. (2006) was the ef- 
fect of blending from faint stars within the seeing disc of a bright 
microlensed source (as opposed to the case for which the faint star 
itself is magnified). We find that for our simulations this is a very 
weak effect. In our reference field, a sample of bright events with 
Is < 17 (i.e. source magnitude, not baseline magnitude) contains 
less than 0. 1 per cent that are heavily blended (fs.m < 0.2) and only 
around 1 per cent with moderate blending (/sj„ < 0.5). For these 
bright event with /s < 17, the error in the recovered tE,4p is negligi- 
ble; for our reference field the value of (/E,4p) is underestimated by 
only 2.7 per cent. It is clear that the problems caused by blending of 
bright source stars will be dwarfed by the other, much more preva- 
lent, cause of blending, i.e. microlensing of faint blended source 
stars below the magnitude cut-off. 



6 THE UNCERTAINTY IN THE OPTICAL DEPTH 
MEASUREMENT DUE TO BLENDING 

It is important to quantify whether the assumption that bright mi- 
crolensing events are unblended leads to any bias in the measure- 
ment of the optical depth. As has been shown above, this assump- 
tion can lead to a significant error in the value of the mean event 
timescale. However, it has been argued that this effect cancels out 
due to an overestimate in the number of events with bright source 



stars (e.g. Afonso et al. 2003). It is the misidentification of the true 
baseline source magnitude that causes this overestimate, since a 
sample of bright events will be contaminated by fainter source stars 
that appear brighter than the magnitude cut due to blending. In this 
section we test the validity of this assumption by taking a sample 
of bright microlensing events from our simulations and calculating 
both the input and output optical depth (i.e. both the true value and 
the measured value as calculated using the conventional methods). 

6.1 Method 

6.1.1 Input optical depth 

The first step in this procedure is to calculate the input optical 
depth, Ti„, for bright source stars (i.e. /[„ < 17). The optical depth is 
the probability that a given star is magnified by greater than 3 / VS. 
This can be estimated for a stellar field by counting the total length 
of time that any of the bright stars are magnified by greater than 
3/ Vs, divided by the total number of bright stars and the total du- 
ration of the experiment. Since the length of time for which any 
microlensed star is magnified by greater than 3/ VS is simply given 

by 2fE,in .^1 - "Qi„, the optical depth for our simulated catalogue 
can be calculated using the following equation (e.g. Popowski et 
al. 2005), 
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Figure 7. The ratio of (efficiency corrected) fitted vs input mean event timescales for our simulated fields as a function of stellar density. The thin and thick fines 
correspond to 5-parameter (i.e. unconstrained blending) and 4-parameter (i.e. /s = 1) fits, respectively. The dotted/solid/dashed lines denote low/medium/high 
quality seeing. Note that the 5-parameter fits are very stable and therefore the lines are practically coincident, i.e. all have (?E,oLt)/{fE,in) ~ 1- 



(6) 



where T is the duration of the microlensing experiment, A'j, i„ is the 
number of stars in the input field with /;„ < 17 (i.e. before the field 
has been convolved to ground based seeing), and A'cvcnts is the total 
number of microlensing events that have been generated. Note that 
A'events corresponds to the number of events that have been gener- 
ated, not the number of events that have passed our detection cri- 
teria. Since we are interested in the optical depth to bright stars, 
A^events Corresponds to the number of events with source magni- 
tude brighter than this cut-off, i.e. /[„ < 17. Also, we only consider 
events that have peak time within one of the observing seasons, i.e. 
T does not include the off-season. Therefore our simulations, which 
are based on three year OGLE-Il observations, have T = 778 d. 



6.1.2 Output optical depth 

The input optical depth determined from equation l|6) must then be 
compared to the observed optical depth, Tant, which can be calcu- 
lated using the conventional equation (e.g. Udalski et al. 1994a), 



DoPHOT once the field has been convolved to ground based see- 
ing) and A'events is the total number of microlensing events that have 
passed the event detection criteria and have well-constrained val- 
ues for the fitted timescale (see Sections [3l & 15. 1 .B . A'events "'^'y 
incorporates events with baseline magnitude l^nt < 17 and impact 
parameter i/o.out < 10. The detection efficiency corresponds to the 
fraction of events that pass our event detection criteria (see Section 
I5.1.3t . Note that equation ^ is not the only way one can estimate 
the optical depth; for example, Popowski et al. (2005) advocate a 
slightly different approach that incorporates the maximum amplifi- 
cation into the detection efficiency. 

We estimate two different values of t: the optical depth es- 
timated under the assumption that bright events are unblended 
(Toiit,4p) and the value estimated when blending is incorporated into 
the fit (Tout.sp)- Each of these t estimates are based on samples of 
events created using the respective detection critera, e.g. for rout,4p 
the detection criteria requires Ko4p < 1, while for To„, 5p the detec- 
tion criteria requires wo.Sp < 1- Note that for our blended fits, we 
only retain events for which the source magnitude is brighter than 
^s,5p =17 mag, not the total baseline magnitude. 



f(?E,out,i) 



(7) 



where A'^ om is the number of stars in the output field brighter 
than /(,„, = 17 mag, (i.e. the number of bright stars detected by 



6.2 Results 

To assess the reliability of recovering t for our two methods, we 
calculate the ratio Tou,/Ti„ for each of our simulated fields. This 
information is given in Table[T]and is shown in Fig.[8] 
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In general, the 5-parameter blended fits produce good agree- 
ment between Tom and Ti„, with an error of less than 20 per cent. 
Note that the error is less than 10 per cent for all fields except those 
with lowest quality seeing (i.e. 2. 1 arcsec). 

This behaviour can be understood as follows. The dominant 
effect for the 5-parameter fits is the error in the estimation of 
This value is determined from the number of observed stars with 
^ut < 17 mag, but since a number of stars on the image are blended 
we find that that N-^ ^^t > ^*,in- Blending can affect the recovery of 
A'star in two ways: firstly, faint stars can be brought into the bright 
regime by being blended with other faint stars; and secondly, a pair 
of bright objects can be unresolved and hence be detected as only 
one object. These are clearly two competing effects, yet the fact that 
^*,out > A'^^in shows that the former effect is dominant, reflecting 
the steepness of the luminosity function at the main sequence turn- 
off. In general we find that Af*.out overestimates the number of stars 
by around 10 per cent. However, for the fields with worst seeing 
the error in A'* o^, can reach as high as 36 per cent, which explains 
why the error in Tom increases for fields with worse seeing. In the- 
ory one could correct for this issue using simulations, although in 
practise this may be difficult to estimate accurately due to certain 
observational effects (Sumi et al. 2006). 

However, even if we could correct for this issue the problem 
is not fully resolved. For example, if we use N^-in in equation l[7j 
the 5-parameter fit still overestimates the optical depth by around 
10 per cent for all fields. This is in part due to minor effects that 
influence the recovery of Tom, such as incorrectly estimated source 
magnitude and impact parameter from the model fits. In addition to 
this, there is an error associated with the fact that to obtain equa- 
tion l|7j from equation one assumes that the impact parameter 
is uniformly distributed. In practise this is not entirely true since 
events with impact parameter close to 1 have lower detection effi- 
ciency than events with larger magnification (i.e. smaller z(o). This 
problem can be overcome by using an efficiency that is a function 
of both /e and mq (e.g. Popowski et al. 2005) or by introducing a 
correction factor determined from Monte Carlo simulations (such 
as those presented in this work). 

The case for the 4-parameter fits is somewhat less subtle. As 
has been stated above in Section 15.31 when a sample of bright 
events are erroneously assumed to exhibit no blending, the fitted 
timescales are shifted towards smaller values (see Fig. |5] and Ta- 
blefTJ. However, the effect on the optical depth is counterbalanced 
by the increase in the number of observed events, i.e. a number of 
events with source star magnitudes fainter than the cut-off are in- 
cluded in the sample due to blending (see, for example, Alcock et 
al. 1997b). For example, in our reference field we find that 22 per 
cent of bright events are actually caused by fainter sources with 
^s,in > 17. In our extreme case of lowest quality seeing and highest 
density this fraction rises to 47 per cent. Owing to the cancellation 
of these two effects, from Table [T] we see that the error in the un- 
blended optical depth is less than 7 per cent for all of our simulated 
fields, even for those with highest stellar density or worst quality 
seeing. In most cases the recovered optical depth is slightly lower 
than the input value, but such an error can be considered negligible 
compared to the statistical errors that will be present in any current 
real-life experiment. It should also be noted that there are a number 
of additional factors that affect the final optical depth, such as the 
fact that many events with ;(o < 1 have kq 4p > 1 after fitting with 
a model assuming no blending (see Fig. |5}. As was discussed in 
Section [53] there is the issue of bright source stars being blended 
with nearby faint stars, but we have found this to have only limited 
influence. The problem described above of incorrectly determining 



W^^ oui also affects the 4-parameter fits, although in this case it is not 
the dominant factor. However, this will still bias the optical depth 
towards smaller values. 



7 CENTROID MOTION AND CENTROID-SOURCE 
OFFSET FOR BLENDED EVENTS 



As has been shown in Section [Z2l it is possible to identify blended 
microlensing events through their centroid motion, or equivalently 
the offset between the baseline centroid and the lensed source. 

In Fig. [9] we show the distribution of the offset between the 
baseline centroid and the lensed source for our simulated fields. 



(rs< 



■ rblend)/^bl 
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where Tbicnd and Tsouicc denote the location of the blend and source 
centroid, respectively, and Fbicnd and Fjomce denote the blend and 
source flux, respectively. From Fig.|9]it is clear that for almost all 
fields a significant fraction of events exhibit noticeable offsets of 
more than 0.2 arcsec. Although Ar is reduced for the bright samples 
of stars, there is still greater than 10 per cent with Ar > 0.2 arcsec. 
The amount of offset is particularly severe for the fields with low 
quality seeing; for example, for field sld2, which has a seeing of 
2.1 arcsec, more than 80 per cent of events have Ar > 0.2 arcsec, 
even for the bright sample. 

However, one must interpret this figure with caution because 
one cannot simply deduce that all of the events with large offset are 
heavily blended. For fields with low quality seeing (2. 1 arcsec), as 
many as 50 per cent of bright, unblended (0.7 < fg < 1.3) events 
have Ar > 0.2 arcsec. Fortunately this problem is reduced for fields 
with better quality seeing, such as our reference field (seeing 1.05 
arcsec); of the bright events in this field with Ar > 0.2 arcsec, only 
8 per cent show no significant blending (i.e. /s > 0.5). However, 
the problem remains for low quality seeing, even in the sample of 
bright stars. Therefore, it is apparent that if one wishes to use such 
information to quantify the level of blending, then this must be used 
with caution. In particular, it is unlikely that any meaningful results 
can be obtained for cases of low quality seeing (e.g. 2.1 arcsec). 

It is also possible to calculate the maximum centroid motion 
(Ar') that would be observed for the events in our simulations. 
Given that the location of the centroid can be expressed as. 



r(0 



Tblcnd^bl 



eA(0 



Ml) 



(9) 



where A(f) denotes the magnification at time /. Given this formula it 
is trivial to derive the expression for the maximum centroid motion, 
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where Apeak is the maximum magnification. This equation can 
equivalently be expressed in terms of the centroid-sourse offset (Ar, 
given in equation[8}. 
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^blend ~^ Apeak ^source 

Clearly this maximum centroid motion corresponds to the differ- 
ence between the location of the centroid at baseline and at peak 
magnification, and hence is an upper limit to the observed motion 
since the event may not be observed exactly at peak magnification. 

The cumulative distribution for Ar' is shown in Fig. [10] As 
expected, this figure shows similar behaviour to the correspond- 
ing plot of baseline centroid-source offsets, although the offsets are 
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Figure 8. The ratio between the input and output optical depth for 4-parameter (i.e. /s = 1; thick) and 5-parameter (i.e. unconstrained blending; thin) tits for 
our simulated fields. The dotted/solid/dashed lines denote low/medium/high quality seeing. 



somewhat smaller. For example, for our reference field (medium 
seeing and density) 30 per cent of all events and less than 10 per 
cent of bright events have Ar' > 0.2 arcsec. 



8 DISCUSSION 

This paper has provided a detailed investigation into the nature 
of blending in gravitational microlensing experiments towards the 
Galactic bulge. As we have discussed in Section |2l it is clear that 
for some bright events blending can be significant. Our simulations 
also support this claim. In Section|4]we showed the true underlying 
blending distributions for our simulated fields, all of which exhibit 
blending to varying degrees. We can conclude that the dominant 
source of blending is from lensing of faint source stars, rather than 
lensing of bright source stars blended with nearby fainter stars. 

However, the results from our simulations (Sections [3]-|6f in- 
dicate that this issue may not be as troublesome as one might fear. 
Although the event timescale is unquestionably affected, the opti- 
cal depth determinations in most cases are robust across a range of 
seeing and density. The fact that the optical depth is reliably recov- 
ered even from a sample of blended events is particularly important; 
however, this reliable recovery appears to rely on a coincidental 
cancelling out of two factors (an underestimation in the timescale 
and a corresponding overestimation of the number of events). This 
finding is not new, but to discover that it holds for various values 
of seeing and density in our simulation is reassuring. Clearly we do 
not advocate that our results should be applied directly to any real 



experiments, since the details will depend on the individual exper- 
imental set-up and event detection criteria. 

The results in Section[5]are of slightly greater concern in rela- 
tion to the recovery of the event timescale. If we (erroneously) as- 
sume that our bright events are unblended then this clearly results 
in an underprediction of the event timescale of between 10 and 40 
per cent. It has been argued that the extent of this problem can be 
reduced if one only deals with high amplification events. For ex- 
ample, the recent papers of Popowski et al. (2005) and Hamadache 
et al. (2006) advocate only using events with amplification greater 
than 1.5 and 1.6, respectively. If we apply this restriction to our 
simulations then the level of discrepancy between the fitted and 
input timescale is reduced slightly; for example, if we apply the re- 
striction that the A > 1.6 for our reference field we find that the 
ratio between the fitted and input becomes 0.81, which is only 
a slight improvement compared to the ratio 0.77 that was found 
for events with A > 1. Although this does improve the situation, it 
is far from resolving the problem. This underprediction could help 
to explain the discrepancy in between Popowski et al. (2005; 

= 15±15d)andSumietal. (2006; (fs) = 28.1±4.3 d), since the 
former work assumes that their bright events are unblended while 
the latter work incorporates blending into their model fitting. How- 
ever, there are large errors on {t-g) for Popowski et al. (2005) and 
hence the discrepancy is only very weak. It is also interesting to 
note that another recent clump giant survey that assumes no blend- 
ing (Hamadache et al. 2006; {t-g) = 28.3 ±2.8 d) finds a value that is 
in good agreement with Sumi et al. (2006), although owing to their 
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Figure 9. The distributions for tlie offset between the baseline centroid and the lensed source for microlensing events in our simulated fields. The upper panels 
show the distiibution for a bright subset of events (laut < 17), such as red clump giants, and the lower panels show the distribution for all events (loat < 19). 
The left panels show fields with differing stellar density (dotted/solid/dashed correspond to high/medium/low density) while the right panel shows fields with 
differing values for the seeing (dotted/solid/dashed corresponds to low/medium/high quality seeing). Note that in our simulations one pixel corresponds to 
0.417 arcsec. 



larger spatial coverage one might expect them to find a larger (fg) 
(e.g., Wood & Mao 2005). 

An aspect of blending that has been ignored in this work is that 
of colour. Clearly if the colour of the blend and the source are dif- 
ferent, then the microlensing amplification will be chromatic (e.g., 
Kamionkowski 1995; Buchalter, Kamionkowski, & Rich 1996). 
The strategy for the OGLE experiment has been to take practically 
all observations in the /-band, with only very limited sampling in 
the V-band for determining baseline colour information. However, 
EROS and MACHO both have colour information, although in their 
optical depth papers (Popowski et al. 2005, Hamadache et al. 2006) 
colour information was not used to discriminate individual blended 
events (although the former work did note that the colour variation 
of their sample was consistent with expectations for a sample of 
unblended events). 

Although the work presented here has dealt with the issues re- 
lating to blending from stars that are coincidentally located close 
to a microlensed source, there is an additional factor that we have 
not considered, namely the luminosity of the lens. By definition ev- 
ery microlensed source must be coincident with a lens star. Various 
works have attempted to address how this fact affects the amount of 
blending and the recovery of parameters from a theoretical stand- 
point (e.g., Kamionkowski 1995; Nemiroff 1997; Han 1998). In ad- 
dition, for one event towards the LMC (Alcock et al. 2001b) and 
at least one event towards the Galactic bulge (Kozlowski et al., in 



preparation) the lens and source have actually been resolved, im- 
plying that for these events the lens contributes a finite amount of 
flux to the total blended light. If more events can be found for which 
the lens and source are resolvable, it will provide a promising way 
to empirically quantify the effect of lens blending. Since the Monte 
Carlo analysis presented here (Sections[3]-[6} has not incorporated 
this lens flux, our work can be considered to underestimate the ef- 
fect of blending. 

One obvious deficiency in this work (and all other similar 
work) is the lack of available deep luminosity functions for the dif- 
ferent lines of sight towards the bulge. Deep HST luminosity func- 
tions, such as the Holtzman et al. (1998) one, would be very valu- 
able if they were available across different bulge fields. A project to 
tackle this issue is currently underway (Kozlowski et al., in prepa- 
ration). 

Data from projects such as OGLE-II can (and have) been used 
for additional purposes not directly related to microlensing, such 
as producing extinction maps (e.g. Sumi 2004) and constraining 
the parameters of the Galactic bar through red clump giant num- 
ber counts (e.g. Rattenbury et al. 2007). These studies will also 
be affected by the problems of stellar crowding. For example, any 
work that wishes to determine the centre of the red clump giant re- 
gion will find it systematically shifted towards brighter magnitudes 
due to blending. From our simulations we can estimate this effect; 
by fitting our observed luminosity function with a power-law plus 
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Figure 10. As Fig.|9]excepl Ihis shows the centroid motion for microlensing events in our simulated fields. The centroid motion corresponds to the motion that 
would be observed if the event were observed at the peak magnification (see Section|7}. 



Gaussian (where the Gaussian represents the red clump stars; see, 
for example, equation 4 of Sumi 2004) one can easily determine the 
shift. For our reference field we find that the Gaussian is centred on 
an I-band magnitude of 16.062 ± 0.004 and 15.918 ± 0.005 for the 
input (i.e. HST) and output (i.e. after convolving to ground-based 
seeing) luminosity functions, respectively, i.e. a shift of 0.144 mag. 
Such a shift can have an important effect in applications such as 
those mentioned above. However, for the issue of extinction, one is 
often more interested in the colour of the centre of the red clump 
region. Since we do not have any colour information in our simu- 
lations we cannot address this point, but from Fig.|2]it can be seen 
that if the blend is caused by a star that is less than ~ 2 magnitudes 
fainter than the red clump region, then there will be a systematic 
blueward shift. 
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In conclusion, we confirm that blending has only a limited 
impact on the recovery of optical depth for the current generation 
of microlensing experiments. In Section|6]we have shown that the 
recovered optical depth is probably reliable to within ~ 10 per cent. 
However, it is now becoming clear that with the improved accuracy 
of future experiments, microlensing surveys will not be able to hide 
from the issue of blending. Such work must undertake a detailed 
and thorough treatment of this effect otherwise their results could 
be subject to significant bias. 
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Table 1. Quantitative comparison of all nine simulated fields. The first four columns describe the setup of our simulated fields, while the remaining columns show the fraction of blended events (columns 6-9), the 
ratio between the input and fitted Je (columns 10 & 11) and the ratio between the input and fitted t (columns 12 & 13). The subscript 4p indicates a four-parameter fit (i.e. /s = 1). while the subscript 5p indicates a 
five-parameter fit with unconstrained blending (see Section |5.1.2) . The fractions of blended events correspond to the input blending distributions, not the fitted distributions. Note that /s > 1 corresponds to negatively 
blended events. Stellar densities are quoted for stars with / < 17. 
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Table 2. Percentage of blended events as a function of /o„t for our field with 
medium quality seeing and medium density. These percentages correspond 
to the input blending distributions (/s.in). not the fitted distributions. 



